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^ ' We apply the GLV Reaction Operator formalism to compute the effects of multiple 

QQ ' elastic scatterings of jets propagating through dense matter. We derive the elastic Reaction 

, Operator and demonstrate that the recursion relations have a closed form solution that 

' reduces to the familiar Glauber form. We also investigate the accuracy of the Gaussian 
dipole approximation for jet transverse momentum broadening. 

o ■ 

^ , PACS numbers: 12.38.Mh; 24.85. +p; 25.75.-q 



O ; I. INTRODUCTION 



New experiments at the Relativistic Heavy Ion Collider (RHIC) will provide high statistics high pr mea- 
surements of photon, lepton and hadron production in A + A reactions [ 0, 0] . These observables provide 



> 

k> , novel tests of pQCD multiple scattering theory. Quantitative comparisons [ H:^J5| to data require under- 

j_j ■ standing of nuclear effects such as parton shadowing [ |j , the Cronin effect [ 0^ ^ ||, p3 a nd non-abelian 

jet energy loss [ |l^, ^ In this paper we apply the Reaction Operator method [ |13|, developed to 

solve the radiative energy loss problem, to derive the well known multiple elastic Glauber distribution [ 14 
in order to illustrate the power of this recursion technique. The summation to all orders in the opacity, 
X — J cr{z)p(z)dz = L/X, of the matter is given in closed form. In addition we evaluate the accuracy of the 
dipole approximation, leading to the Gaussian approximation to Moliere scattering, for the case of moderate 
opacities x < 10 relevant for applications in nuclear physics. 

We describe multiple jet scattering as a series expansion in x, the mean number of scatterings that a hard 
parton undergoes along its trajectory [ Jet production is modeled here by an initial wave packet 

j{p) for a spinless parton in color representation R prepared at finite time to and centered at Xq — (zo,xo). 
The momentum space amplitude in the absence of interactions is then given by 

Mo = *e'f-«j(p)ld..xdH , (1) 

where the dfj dimensional unit matrix accounts for the jet colors. Multiplying |AfoP by the invariant one 
particle phase space element (i^p/((27r)'^2|p|) and taking the color trace one arrives at the unperturbed 
inclusive distribution of jets in the wave packet 
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In the presence of nuclear matter the multiple interactions of a fast projectile propagating through the 
medium are modeled by a scattering potential assumed to be of the form 

K = i;(g„)e^«""" (R) ® r,„ (n) - 2^(5(<70)i;(q„)e-^^"-^" r,„ {R) ® (n) , (3) 

where T^^ (i?) is the generator of SU{Nc) in the dimensional representation of the jet, and T^^ (n) is the 
corresponding generator in the dn dimensional representation of the target. The scattering center is localized 
at position x„, and f (qn) is the Fourier transform of the spatial part of the potential. We here use p to 
denote the 3D spatial part of a vector p and p represents its 2D transverse component. While there are no a 
priori limitations on the functional form of the scattering potential and the corresponding differential elastic 
scattering cross section 

daei{R,T) ^ CrC^{T) |«(q)p 

(i2q dA (2^)2 ' ^ ' 

we here consider potentials with a range ^ A, the mean free path of the jet in the medium. This 
condition together with the hard jet approximation ensures the path ordering of the successive scatterings 
along the jet trajectory. Diagrams with crossed momentum transfers (at points typically separated by 
Zn — Zn-1 — A) are suppressed by factors ~ exp(— A/^). One particular example used in Refs. is the 

Yukawa color-screened potential 

v{(U) = ^y- — 2 = 7 r~ — T7 ■ — A ' 

q« + ^J■ (Qnz + IfJ-n ) [qnz - IfJ-n ) 

where /Lt^ = /i^^ = + Qn a-^d A^ > 1 . 

This paper is organized as follows: in Sec. II we present a systematic recursive way of keeping track of the 
multiple collision Feynman diagrams in terms of "direct" and "virtual" interactions. We discuss the color 
and kinematic structure of the single Born (direct) operator D and the double Born (virtual) operator V. In 
Sec. Ill the corresponding elastic Reaction Operator R ~ D^D + V'^ + V is computed. We find a recursion 
relation for the observed inclusive jet distribution at any fixed order in opacity and use impact parameter 
space resummation to recover a classical Glauber-like formula. In Sec. IV we discuss the implications of our 
results in relation to the frequently employed Gaussian approximation of transverse momentum distribution 
of jets. We show, however, that for x < 10 the Gaussian approximation fails to reproduce the power law 
tails of the jet pr distribution as well as the ~ logQ^^,^, enhancement of (p^)- 



II. TENSORIAL BOOK-KEEPING AND DIAGRAMMATIC CALCULUS 



The interaction Hamiltonian in a medium with scattering centers in the high jet energy spinless ap- 
proximation is given by 

N 

Hj{t)= / d3s^t-(x-xi)T,(i)0t(jJ,t)T,(i?)i5(t)0(x,i) , (6) 

where D{t) = i dt- Each interaction of the colored parton propagating through the nuclear matter results 
in one power of the scattering potential Eq. (^ in the diagrammatic expansion. There is no limit on the 
number of scatterings the jet can have at a particular scattering center at position z„. Using the terminology 
of Ref. [ |l^ we call such interactions single Born, double Born, e.t.c. It can be seen from Eq. (||) that two 
momentum transfers resulting from two interactions with the scattering potential are needed to produce 
one power of the elastic scattering cross section a^i ■ In performing the average over the transverse (relative 
to the jet trajectory) position of the scattering center (• • •)^ = J d'^h/A±{ • • • ) one order in the optical 
thickness (opacity) Naei/A± = L/X is thus generated in three ways: a.) single Born or "direct" interaction 
in the amplitude and the complex conjugate amplitude; b.) double Born or "virtual" interaction in the 
amplitude and no interaction in the complex conjugate amplitude; c.) no interaction in the amplitude and 
double Born interaction in the complementary complex conjugate amplitude. More momentum transfers at 
a fixed longitudinal position Zn, e.g. a triple Born interaction, will produce 0(0^) corrections to the opacity 
expansion and are here neglected. 
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It is therefore sufficient to consider single and double Born terms in the opacity expansion approach, i.e. 
when the jet propagates by a scattering center at position Zn it can either miss the center, interact once, or 
interact twice as illustrated in Fig. 1. 




p,c 



A,- 



ii,...,i„.i 



p qn.an 



FIG 1. The diagrammatic representation of the interactions with a fixed scattering center located at z„ 
up to the double Born term ( /„ Aii-.-i^_-^^, Dn -4ii...i^_j, and Vn Ai;^.-.i^_-i ) are shown. 

The kinematic and color structure of the scatterings is contained in the unit ( / ) , direct ( D ) and virtual 
{V ) operators. The path ordering of the projectile interactions in the medium in the high energy eikonal 
regime facilitates the introduction of a convenient tensorial book-keeping notation that classifies the diagrams 
according to the type of scatterings that the parton has experienced along its trajectory. At order n in the 
opacity expansion there are 3" possible amplitudes classified by a set of n indices ii - ■ - in- Here im = 



indicates that the jet misses the m target 



1 means that it scatters once, and i„ 



2 indicates that 



it scatters twice. The diagrams in Fig. 1 can be algebraically represented as follows: 



■^ii ■■■'in 


-i,o(p, c) 






■^ii ■■■'in 






-lb: C) 


-^ii ■■■'in 


-l,2(P, C) 


— ^-^ii ■ ■ -in 





(7) 

The set of amplitude indices in Eq. (|^) thus encodes the complete history of the jet interactions. Repeating 
the basic operator steps in Eq. (|^) any amplitude that includes parton scatterings inside the medium can be 
iteratively derived from the unperturbed jet production amplitude 



m— 1 



+ Dm + 62,i^ Vm Mo (p, c) 



(8) 



Time (or longitudinal coordinate) ordering is implicit in the above formula for the high energy eikonal limit 
under consideration here. 

To proceed we need to know the color an kinematic structure of the direct and virtual interactions. We 
use the notation A{p) = (p^ and T{p) = for the momentum space representation of the propagator 

and the vertex factor. The vertex factor is treated here in the high energy spinless limit. The direct iteration 
step (with Zn > Zri-i > ■••) in Eq. (^ reads 



A,...,^_,ip-qn,c)A{p-q„)viqn)r{2p~qn) e"^"'^^"-''"^ TaMTajR) , (9) 
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where we have taken into account that the amphtude of order n—l has to be evaluated at momentum p — q„ 
if the jets emerges on shell with momentum p. The double Born amplitude at the same external potential 
is similarly given by 

Ai--j„_i,2(p, = y (2^ -^n-'i^-i ip-Qn- q'n, c)A{p - g„ - qn)T{2p - 2ql^ - x 

X A{p - g;)r(2p - q'Jviq'„) e*(9"+«:.)("o) (^)r^^^ (^) j^^^ (^) . (iq) 



III. ELASTIC RECURSION OPERATOR 



The observed jet inclusive distribution dN{p,c) = X]^o '^-^'"HPi c) is expanded in the opacity series, 
where the contribution to order is given by 



dN^"\p,c) = A^^■■■^"{P,c)A^■■■^AP,c) ^C^TtJ2---T. A^'-'"iP,c)A,...,AP,c) • (11) 

ii=0 i„=0 

The amplitudes A'^^"''^"{p,c) in Eq. ( pi] ) are the complementary amplitudes given by 

n 

^^--"(p,c) = Mt(p,c) n [So,^yi + h.^Dl + <52,,„] . (12) 

m— 1 

The trace is taken over the color matrices and the binomial coefficient = Nl/nl{N — n)! w N'"-/nl is 
introduced to take into account the number of combinations of n scattering sites out of the total N. 

Performing the sum over the first n—l interaction points and using Eqs. (|8|,p^ we obtain a simple recursion 
identity which relates dA^'"' to dA^'"^^^ through the Reaction Operator R 

d7V(«) ^ ^*--'"-ii?„Ai....„_i , Rn = DiDn + Vn + Vl ■ (13) 

Consider first the direct part dA^^"-* (Dir.) oc A!'^"'^'^-^t)'^I)Ai^...i^_^ of the Reaction Operator. Performing 
the color traces we take into account the identity Tr Ta(i)T(,(j) — 5ij5abC2{T)dT / dA, which enforces that the 
positions of the ordered scattering centers in the ainplitude and its complementary are identical. We use the 
form of the scattering potential specified by Eq. (^ as well as the the high energy eikonal approximation 
where the deviations due to the in-medium interactions from the initial jet trajectory are small and £'+ ~ 2E. 

dA^(")(Dir.).C- j %"^)3^"' X 



E^flzn -ql^-cg^ + ie E+q'^^ - q'^\ - q;2 - ie 

The qzmq'zn integral can be performed by closing the contour in the lower/upper half plane. We note 
that due to the shor t range of the sc attering p otential relative to the mean free path the residues from the 
qzn — ~i\/^n + ^-iid q'zn — +i-\/q^^ + poles are exponentially suppressed and only the residues at 
qzn ~ — «e + ^n/E'^ and « +ie + q'^^/E'^ contribute. In the high energy limit i?+ » ^ > ^^JE^ the 
quadratic terms from the residues in the exponent can be neglected. This removes the correlations between 
the scattering centers and there are no coherence effects at the jet level. We note that this is an important 
difference from the case of induced gluon radiation discussed in Ref. [ |l3| . Emitted gluons are much softer 
than the parent parton and non-negligible (j) ^ (x.,- — xq) ■ (k— Q)2/fc"'" phases control the non-abelian analog 
of the Landau-Pomeranchuk-Migdal effect [ M. Taking into account that the Glauber thickness function 
(at a fixed impact parameter bo) ^(bo) = jdz p{ho, z) — N/A± is expected to vary slowly with impact 
parameter: 

^-.(.„-0-(b-bo)^^ -^(2^)^^^(q„-q;) . (15) 
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Eq. (^_5[) is a key simplifying relation valid in the case \fA\_ ^ 1/ yL because it diagonalizes Eq. (|lj) in the 
transverse momentum variables. 

This leads to a simple recursion relation for the Direct contribution to the jet spectrum to the momentum 
shifted distribution of order n — 1 . 



ciiV(")(Dir.) -C,^ / d2q„^*i 



dcTel {R, T) T(bo) ,q„.bt ,q„.b 



Ai - i„_i(p,c) , (16) 



where h — i Vp is the impact parameter operator conjugate to the transverse momentum p acting to the 
right, and b''' = —i Vp its Hermitian conjugate acting to the left. 

Next consider the virtual contribution dA^(")(Vir.) = A'^'""^"-^ {V + V^)Ai^...i,^_^ of the Reaction 
Operator, Eq. (|l|). 



d7V(")(Vir.) = 2Re 



A' " nwAi...j„_i(p- gn - g„) w(qn)w(q„) -. X 



(27r)3 



(27r)3 



dA 



,-i(q„+q;j-(x„-Xo)\ ^ ^^^^ 

One notes that the momentum space integrations in this case could have been performed already at the 
amplitude level, Eq. (|l0|). Using the same set of approximations as for the Direct part we first perform 
the qz n integral as well as the impact parameter average and note that the latter constrains q„ + q^^ = 0. 
Picking the residue at qzn ~ ^Qzn ^ 

dq'z^d"^ qn 



dAr(")(Vir.) = 2Re 



{2n) 



A' 



T{ho) CrC2{T) 



N 



dA 



E+q'.n 



„/2 

n 



q;2 + ie 



(18) 



It was shown in Ref. [ ^ that for a broad class of real spherically symmetric potentials of short range 
<ti A the remaining integral can be performed. The screened Yukawa potential, Eq. (|3|), belongs in this 
category when fiX ^ 1. We note that the contour integration produces a real result and a factor of ^ that 
cancels the factor of 2 arising from the Virtual terms in both the amplitude and its complementary. The 
recursion relation for the virtual part of the Reaction Operator can then be written as 



dAf(")(Vir.) = J d\nA''- 



daMT) r(bo) 



(-1) 



d2q„ N 

Combining Eqs. (E6^,n9h one recovers the full structure of the elastic Reaction Operator 



(19) 



dN(-\p,c)=C^ 
1 



da,iiR,T) T{ho) 



d^qn 



-icin-h' 



N 



»q„ b _ ^ 



•^ii---i„_i {P: c) 



T(bo 



daei {R, T) 



-q, Vp _ I 



diV(0)(p,c) , 



(20) 



where we used the /N"^ « 1/n! (Poisson) approximation. 

Eq. ( po| ) is the main result of this paper and has transparent physical interpretation illustrated in Fig. 2. 
The first term in the kinematic structure of the Reaction Operator corresponds to the direct elastic scattering 
part. The next two terms have delta function strength in the forward direction and correspond to the virtual 
corrections (in the amplitude and its complementary complex conjugate). These terms enforce unitarity in 
the GLV Reaction Operator formalism. Given any initial parton distribution one can recursively solve from 
Eq. ( pO| ) for the final inclusive distribution of jet that have penetrated the medium characterized here by 
its Glauber thickness function T(bo). The method used here is quite general and has several important 



5 



advantages. Within the eikonal approximation it provides exact answers to problems related to multiple 
elastic and inelastic processes in QCD media. In particular it was applied to the case of gluon radiation^ 
from jets produced in heavy ion reactions [ For more details on its uses in jet tomography see Refs. [ 
^, P . Since the solution is known to all orders in opacity, cases ranging from thin media with few scatterings 
to thick media with formally very large number of scatterings can be studied. 



R = 









■i ' 








FIG 2. The graphical structure of the Reaction Operator R = & D + V + V^ is illustrated. It represents 
all possible {t = co) on-shell cuts through a new double Born insertion. 



In the case of multiple elastic scattering, the reaction operator is so simple that it is possible to sum over 
all orders in closed form. As usual it is most convenient to perform the summation in impact parameter 
space. The distribution in impact parameter space conjugate to the transverse momentum p is 



oo -. 

dN{h) = V i- [r(bo) (ae;(b) - a,i)r dN(°\h) = e^(^«) (*-(^)---')(iiV(«) (b) 

^ — ^ n ' 



(21) 



n=0 



where the Fourier transform of the differential cross section was denoted by {2Tr)'^ da ei/d'^q{h) = (Tei(b) 
and o'e;(0) = cr*°*. The impact parameter space representation Eq. ( pT|) depends only on the difference 
aei(b) — deiiO). Transforming back to momentum space, we recover the parton version of the Glauber 
multiple collision series [ n 



dN{p) 



n=0 



f[d\.-p^dN(^\p-ci, 
^J- a el d^m 



(22) 



where x = o'e;T(bo) = L/X. Eq. (22) has the familiar physical interpretation of a Poisson random walk in 
transverse momentum space distributed according to da^ I o^i ■ 



*In the case of medium induced gluon bremsstrahlung the Reaction Operator R retains its general operator form 
Eq. ( |l^ ) but its color and kinematic structures are different since the on-shell t = oo cuts (see Fig. 2) go through 
both the jet and gluon lines that interact with the double Born term. 
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IV. DISCUSSION 



Impact parameter space resummation is often used in conjunction with small impact parameter approxi- 
mation of the differential cross section, Eq. (^). We show below that although such an approach is analytically 
appealing, it gives jet distributions that may differ substantially from the exact formula. 

The Fourier transform of the normalized cross section is given by: 



where b = |b| and in the quadratic term in Eq. ( |23| ) the log 2/(1.08 /i 6) multiplicative factor has been 
absorbed into a 6-independent constant ^ due to its small logarithmic variation. We consider the transverse 
momentum broadening [ ^ |l^ of a jet propagating in the direction, i.e. dN'^^'^ /(Pp = <^^(p)- In the 
approximation that dominantly small impact parameter scatterings are important, the momentum space 
distributions from Eqs. ( ^ , ^2|) reduce to the classic Moliere form [ ]l5| 



rf^be'P''-— ^ ^ . (24) 

(2^)2 xf^H 27T xi^H ^ ^ 

The resulting distribution is of Gaussian form and has a width of x i-e. (p^) = X 




p,[GeV] p,[GeV] 

FIG 3. The final jet pr distribution is shown versus pr for two diflferent opacities x — 3 (Fig. 3a) and 
X = 10 (Fig. 3b). We compare the full result (without the delta function contribution at pr ~ 0) to the 
Moliere Gaussian approximation with ^ = 1 and ^ = logx- In this example we use fi^ — 0.25 GeV^. 



To asses the accuracy of the Gaussian ansatz for finite x relevant in nuclear physics, we compare the 
Fourier transform of Eq. (^) to Eq. (^. Note that in going from Eq. (^ back to momentum space, a 
g-x S'^{p) component arises. This component is of course not displayed, and for finite p numerical integration 
converges much more rapidly if unity is subtracted from exp[CT(b)r(bo)] in the integrand (|2^). 

In Figs. 3a, 3b numerical results are shown for opacities x = 3 and x = 10- For large opacities x ^ 10 
(see Fig. 3b) the Gaussian approximation gives better results than for small ones x ^ 5. However, as is well 
known, at larger transverse momenta it fails to account for the power law pr tail of the inclusive distribution. 
Unlike familiar Moliere scattering in atomic matter, there are no form factors to limit the growth of the high 
Pt Rutherford tail out to the kinematic limit Q^ax -^oM in the case of QCD. Our results for the rms 
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transverse momentum kick shown in Fig. 4 demonstrate that in semi-hard and hard pQCD processes there 
is a logarithmic in Q enhancement of the mean squared transverse momentum, (p^) ~ x log(l + Q^/xM^) 
that is missed in the dipole approximation. Of course for large Q radiative energy loss also contributes 
substantially to the pt broadening and should be taken into account in applications. 
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FIG 4. The transverse momentum broadening of jets is shown versus Qmax for two different opacities 
X = 3, 10 as in Figs. 3a, 3b. We here use fj,^ = 0.25 GeV^ for illustration. 



In summary, we showed that the Reaction Operator formalism developed in [ easily recovers the well 
known Glauber limit for elastic multiple collisions in pQCD matter and sheds light on the accuracy of the 
dipole approximation. In the future it would be interesting to combine both elastic and the LPM suppressed 
inelastic processes in this method and apply the results to jet acoplanarity observables suggested in Refs. [ 
SB- 
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